Pdf s of Derivatives and Increments for Decaying Burgers Turbulence 
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A Lagrangian method is used to show that the power-law with a -7/2 exponent in the negative 
tail of the pdf of the velocity gradient and of velocity increments, predicted by E, Khanin, Mazel 
and Sinai (1997 Phys. Rev. Lett. 78, 1904) for forced Burgers turbulence, is also present in the 
unforced case. The theory is extended to the second-order space derivative whose pdf has power-law 
tails with exponent -2 at both large positive and negative values and to the time derivatives. Pdf's 
of space and time derivatives have the same (asymptotic) functional forms. This is interpreted in 
terms of a "random Taylor hypothesis" . 

PACS number(s) : 47.27.Gs, 02.50.Ey, 05.60.-k, 

I. INTRODUCTION 

E, Khanin, Mazel and Sinai pt] made various predie- 
tions eoneerning the one-dimensional Burgers equation 



IV, we derive 



dtu + udxU ~ vd^u + /, 



(1) 



with viscosity v and a random white-in-time force /(x, t) 
which is homogeneous, periodic and smooth in the space 
variable. One prediction concerns the probability den- 
sity function (pdf) of the velocity gradient ^ = dxU. Ac- 
cording to Ref. ||l|, in the limit u ^ Q, the statistically 
stationary solution of (||) has a pdf 



P(0 « 1^1 



-7/2 



for^ 



(2) 



This power-law range is due to preshocks, nascent shocks 
with a cubic root structure, as discussed by Fournier and 
Frisch 0. There has been an interesting controversy 
about this negative tail of the pdf, which we shall not 
try to summarize here (see, e.g., Refs. 0-01 and references 
therein). There is no complete proof at this moment of 
the validity of the -7/2 law, but significant progress has 
been made recently B. We shall not dwell now on the 
issue of the validity of the -7/2 law for forced hurgulence 
(Burgers turbulence). 

It is our intention here to show that the -7/2 law is also 
present in unforced decaying burgulence. Specifically, we 
shall consider solutions of (|l|) in the limit u -^ with 
/ = and random zero-mean-value initial conditions 
uq{x) which are periodic (a unit period is assumed for 
convenience), statistically homogeneous and sufficiently 
smooth. An instance is to take Gaussian initial condi- 
tions with a spectrum decreasing exponentially at high 
wavenumbers. Such "large-scale" initial conditions will 
develop nonsmooth features (preshocks and shocks) after 
some (random) time. 

The paper is organized as follows. In Sec. |ll| we con- 
sider the deterministic problem in Lagrangian coordi- 
nates and identify the preshock events leading to large 
negative gradients. In Sec. HI we derive the -7/2 law 



for the pdf of the first derivative. In Sec. 

similar laws for higher-order space derivatives and the 

time derivative. In Sec. |^ we derive the corresponding 

results for the pdf of space increments. 

make concluding remarks. 



In Sec. VI 



we 



II. THE LAGRANGIAN REPRESENTATION 
AND PRESHOCKS 

In the absence of force and of viscous dissipation and 
as long as no shock has appeared, the Burgers equation 
(0) has the obvious solution 



u{x, t) — UQ{a), a : 



Lj X, 



where 



tuo{a) 



(3) 



(4) 



is called the naive Lagrangian map. This is indeed just a 
statement that the velocity of a fluid particle is conserved 
in Lagrangian coordinates. (Following standard tradition, 
we denote Lagrangian initial coordinates by a and Eule- 
rian coordinates by x.) 

A remarkable property of the unforced Burgers equa- 
tion in the limit of zero viscosity, which follows from the 
Hopf-Cole solution (see Refs. 0j§l f^or details), is that (|^) 
remains valid in the presence of shocks provided the naive 
Lagrangian map is replaced by the (proper) Lagrangian 
map Cf. The latter is defined as follows. First, we define 
the initial potential (up to an additive constant) by 



uo{a) = -dai'oia). 
We then define the Lagrangian potential by 



(p{a,t) 



+ tVo(a), 



(5) 



(6) 



and observe that the naive Lagrangian nrap is simply the 
negative gradient of the Lagrangian potential : 



Lta = ~—ip{a,t). 
The Lagrangian map is defined as 



d 



Cta = ~—ipcia,t), 



(7) 



(8) 



where ipc{a,t) is the convex hull with respect to a of 
the Lagrangian potential ip{a,t). The convex hull of a 
function /(a) can be defined as the smallest piecewise 
differentiable function which is greater or equal to /(a) 
for all a and such that its derivative is nonincreasing. 

The graph of the convex hull oiip{a, t) is made of pieces 
of the graph of the function ip{a, t) joined by linear seg- 
ments, sitting over the Lagrangian shock intervals, as 
shown in Fig. |l|. 



(p(a,t) 




a" a+ \ a 

FIG. 1. Lagrangian potential and its convex hull in the 
presence of a shock interval extending from a~ to a^. 

Hence, the Lagrangian map coincides with the naive 
Lagrangian map except over the Lagrangian shock inter- 
vals where it is constant (see Fig. ||). Thus, (||) with Lt 
given by (0) remains valid outside the Lagrangian shock 
intervals. 



t>t, 





FIG. 2. Lagrangian map before (a) and after (b) the ap- 
pearance of a shock. The naive Lagrangian map is shown as 
a dashed line. 

We can now use this solution (||)-(§) to calculate the 
Eulerian velocity gradient, i.e. its first-order space deriva- 
tive, in Lagrangian coordinates. Differentiating (j^) and 
using (Q), we obtain 



dxu{x,t) = Uo(a) 



u'o{a) 



1 



daX I + iuQ(a)' 



(9) 



where u'Q{a) = duo{a)/da. We immediately observe that, 
for t > 0, the only way in which this gradient can become 
large and negative is to have a very small denominator 
in (0). For the kind of smooth initial conditions con- 
sidered here, the denominator is necessarily positive for 
sufficiently small times. Let a* be the location where the 
initial velocity gradient Wg (a) achieves its minimum over 
the period. At 



t = t^, = min 



1 



u'oia) 



(10) 



we have the first preshock, i.e. a shock is born M. Subse- 
quently, other (less negative) local minima of Ugia) may 
also produce preshocks, provided the corresponding lo- 
cation has not already been captured by a previously 
generated "mature shock" . Thus, large negative Eulerian 
gradients must come from the neighborhood of preshocks. 
We now recall, for use in later sections, the local (nor- 
mal) form of the Lagrangian and Eulerian solutions near 
a preshock. Let a* be a local negative minimum of UQ(a). 
We then have (generically) 



Mo(a*) < 0, Uo("*) = 0> Uf/i'^*) > 0- 



(11) 



Taylor expanding near a,, we have 
uo{a) ~ Mo(a*) + Uo(a*)(a — a*) + 



6 



(a — a^)^ 



(12) 

By (y), for t near i* — — l/ug(a*), the naive Lagrangian 
map is given by 

a; ~ a, -I- iuo(a*) -|- — (a — a*) 



6 



(a — a,)"^ 



(13) 



Hence, for given x and t, near a;(a*,i*) and i*, respec- 
tively, the naive Lagrangian map can be inverted by solv- 
ing (to leading order) a cubic equation. For t < t^, 
this equation has a single real solution and the naive La- 
grangian map coincides with the Lagrangian map. For 
t — t^, the time of the preshock, the equation simplifies 
and its solution reads 



6 



-,1/3 



i*u[,"(a*) 



{x — x(a,,i*)) 



(14) 



Substitution in (02^ gives 



u{x, i*) ~ uo(a*) — 



6 



t*UQ"(a,) 



{x — a;(a*,i*)) 



1/3 



(15) 



This is the well known Eulerian cubic root structure of 
preshocks. For t slightly in excess of i», the naive La- 
grangian map is not monotonic and cannot be inverted 



(otherwise there would be three branches) . There is now 
a shock. The corresponding shock interval can be de- 
termined by using the convex hull construction on the 
Lagrangian potential. To leading oder, it is found that 
the shock interval extends from a~ to a"*" which are such 
that x{a^, i) = a* + tuo{a*), namely 



± 



6{t-U) 



-,1/2 



tlu'^'ia,) 



(16) 



This condition expresses that, to leading order, the Eulc- 
rian location of the shock remains fixed in a frame moving 
with the velocity uo(a*). 

In Fig. 0, we have sketched the Eulcrian structure of 
the solution at three times, just before, at, and just after 
the time t, of the preshock. 



u(x,t) 




t<t 



t = t 



t>t. 



FIG. 3. Eulerian structure of the solution (a) just before 
a preshock; (b) at the time of a preshock; (c) just after a 
preshock. 



III. THE PDF OF THE VELOCITY GRADIENT 

Our purpose is to derive the behavior for ^ —^ — cxd of 
the pdf of the velocity gradient 

where angular brackets denote ensemble averaging over 
the random initial condition uq. By the assumed homo- 
geneity, p(^; X, t) is obviously independent of x (and will 
subsequently be denoted p(^,t)). It follows that 



p{tt) = { / si^~dM^,t))dx 



(18) 



Having thus a representation of our pdf as a space in- 
tegral over the Eulerian coordinate x, we can make the 
change of variable from Eulerian to Lagrangian coordi- 
nates, using the map Ct- The same idea was used in 
Ref. m to calculate the Fourier transform of the Eule- 
rian solution during the early phase of regularity (before 
the appearance of shocks). This idea also works for later 
times provided we use the Lagrangian map, which differs 
from the naive Lagrangian map (jj) only by the exclusion 
from the basic periodicity interval [0, 1[ of the Lagrangian 



shock intervals. Let us denote by DL{t) the set of so- 
called regular points, i.e. Lagrangian points which do not 
belong to shock intervals. Using the Lagrangian repre- 
sentation (^ of the velocity gradient, we obtain from ( [l^ ) 



Pi^,t) = 



'DlH) 



^u- 



"o(a) 



l + tu'^{a) 



|l + iUo(a)|da 



(19) 



Note that 1 -f- tu'Q{a) is the Jacobian of the Lagrangian 
map. Hereafter we shall use several times the formula 



Sifiy)) = Y. 



1 



l/'(y.)l 



^iy-Vj), 



(20) 



where the yj 's are the zeros of / and / is assumed to be 
sufhciently smooth. 

Let us denote by bk the (discrete) Lagrangian locations 
where the argument of the delta function in (M) vanishes, 
i.e., which are the roots of 



e- 



u'oib) 



l + two(6) 



= 0. 



(21) 



Using (|20|) and evaluating the derivative of the argument 
of the delta function in ( [l9|) at the point where this ar- 
gument vanishes, we can rewrite the pdf of the gradient 
as 



P{tt) = 






1 



Kibk)\ 



DL{t) 



S{a — bk) da 



(22) 



where the integral over the delta function may be viewed 
as shorthand for the indicator function of D^ (t) (equal to 
one if bk G DL{t) and to zero otherwise). Note that the 
r.h.s. of (P^ ) has a |f| "^-dependence for large ^, if we take 
only into account the first factor. Actually, we shall see 
that the presence of UQ{bk) in the denominator gives an 
additional |^|^^^ factor and that realizability conditions 
give a 1^1"^ factor, so that the pdf will be proportional 
to \^\-y'. 

So far, we have not made any expansion. Let us now 
concentrate on the case of large negative ^'s. As ob- 
served in Sec. O, this happens only in the neighborhood 
of preshocks. The latter originate from Lagrangian lo- 
cations at which uo{a) has an inflection point satisfying 
(|ll|). Let a*j be the discrete set of such locations. 

We show now by perturbation theory that, for each 
such point, there are zero or two roots of (Elf). Indeed, 
using the Taylor expansion (f2h, in (pi|), we obtain 



(6- a^jf 



tu'oia^j) 



-—-l-tu'oia^j) 



(23) 



which has either two roots (denoted 6 ) or none, depend- 



ing on the sign of the r.h.s. Defining t^ 



-l/u'^ia^j), 



it is now convenient to distinguish the cases t < ti,j and 
t > t^j, corresponding respectively to before and after 
the preshock. Before, the Lagrangian shock interval near 
a*j is empty; the two conditions that t < t<,j and that 
the r.h.s. of (23) be positive read 



-<-o(«*.)<-7-^- 



(24) 



After i*j, the Lagrangian shock interval is defined by 
(p^). Since shock intervals are excluded from the inte- 
gral (p2[), acceptable solutions must be outside of such 
intervals. This and t > t^,j gives again two conditions, 
namely 



11 // ^ 1 

-i+2^<"°("*^^<-i- 

We observe now that the two conditions 
may be written as a single condition 



(25) 



and m 



1 /^ N 1 1 



(26) 



which, for large negative ^, restricts UQ(a*j) to being near 
— 1/t in a small interval of length — 3/(2t^^). We shall 
denote by ]lps(uo(a*j);i,^) the indicator function equal 
to one if Wg (a*j ) is in this "preshock" realizability interval 
and to zero otherwise. 

The integral over a appearing in (|2^) is over the com- 
plement of shock intervals. The above analysis takes care 
only of those nascent shock intervals which can be cal- 
culated perturbatively. A further condition is that the 
point a^,j should not be within a mature shock interval 
which was created before time t. Because of the convex 
hull construction, this is a global geometrical constraint 
which cannot, in general, be expressed perturbatively. 
We shall denote by II ^j^ (t) (a) the indicator function equal 
to one if a is outside such a "global" shock interval and 
to zero otherwise. 

We now return to (p2f). Since bk — bf, which is close 



to the a^j^s where Uq vanishes, we can use the Taylor ex- 
pansion |wo(6^)| — \b^ — a^j\uQ'{a.^,j), which by ( p3| ) takes 
the same value for b'J and b~ . Hence, the contribution to 
the integral in ( p2|) of the two bk points in the neighbor- 
hood of each point a,j is 21ps (uo(a,j); t,^) lDG{t){o-*j), 
where the products of indicator functions takes care of 
the shock interval exclusion. Using (Eq), we now obtain 

^1/2 

p(C,i)-^^x 



E 



21lps('«o(a*j);i,C)lDG(t)(«*i) 



{K'(«*j) 



-^^-1-tu 



o[0'*j) J 



1/2 



(27) 



Since the sum is over points of vanishing Uq with u'q > 0, 
by use of (BO) this may be rewritten as 



PiU) 



m 



1/2 



HI' 



'\^MM)!!^W) 






1 — tuo(a) 



1/2 



xR{u'^'{aj)lps{u'o{a)-t,0'^Da{t}{a) 



(28) 



where H(-) is the Heaviside function. 

Interchanging now the mean value and the integra- 
tion over a, we observe that, because of homogeneity, 
the integrand does not depend on a. Hence, the in- 
tegration over a can be omitted. Let us denote by 
P3,o{u' ,u" ,u"'\lDc{t) — 1) the joint pdf of the first 
three derivatives of the initial velocity at an arbitrary 
Lagrangian location, knowing that this location is not 
within a mature shock interval at time t. We can then 
write 



P(e,t) 



m 



1/2 



m' 



du'" X 



I^ ^^, {u"'f/h^.^^{u',QM "\tD^(t) = 1) 



-^-1 



tu' 



1/2 



(29) 



In (|29|) , the variable u' is constrained to remain very close 
to —1/t for ^ — > — oo. Assuming that the density pzfi is 
smooth in its u' argument, we can replace the latter by 
— 1/i and carry out the remaining integration over u' , to 
obtain 



p{£„t)^2^t-^D{t)\£,\ 



-7/2 



e 



(30) 



where 



D{t) ^ p du"'{u"'f'^p^.,o i~\,^y I ^Da(t) 



^ 1 



(31) 



This concludes the derivation of the -7/2 law for decay- 
ing burgulence. The time-dependent constant D{t) is ex- 
pressed in terms of a conditional probability which can- 
not be calculated without solving a global geometrical 
random problem 0. A more explicit form is obtained 
for small t : the condition ]l_Dc(f) — 1 may then be omit- 
ted and the integral ( |3l] ) can be calculated, e.g., in the 
Gaussian case. Indeed, for the kind of large-scale initial 
conditions assumed here if, near some point a*, the ini- 
tial velocity gradient Ug achieves a very large negative 
minimum close to —1/t, the other minima will be above 
— 1/t with a probability very close to unity, so that it is 
close to certain that no mature shocks have been formed. 



IV. HIGHER-ORDER SPACE DERIVATIVES 
AND TIME DERIVATIVES 

We begin with the second Eulerian space derivative 
&iu. The method is rather similar to the one used for 



the first derivative. So, we shall avoid repeating details. 
From (|), (|) and (|), it follows that 



dlu = da{dxu) 



daX 



^da 



"o(a) 



l + tu'^^{a) 



<{a) 



l + tu'^{a) [l + tw[)(a)]3' 



(32) 



Denoting by p^ '{£,^ t) the pdf of 9^u, we have, as before 



p('n^,i) 



'DL{t) 



-5 U 



<{a) 



[1 + K(a)]= 



|l + iUo(a)|da 



(33) 



For Gaussian statistics or, more generally, when the prob- 
ability of very large values of Mo(a) is very small, large 
values of 9^m will be due overwhelmingly to small de- 
nominators. That is, they will again originate from the 
neighborhood of preshocks. Near an inflection point a*j 
of the kind considered in Sec. Ill, we have 



dlu 



ig"(a*j)(a- a^j) 



[l + tuQ(a*j) + |TtQ"(a,j)(a — a^^Y 



(34) 



This is an odd function of a — a, _,- which can achieve both 
large positive and negative values. (Hence, we obtain 
power-law tails at both ends.) 

As before, in (p3) we change from a delta function over 
^ to delta functions over those Lagrangian locations h]^ 
where the second space derivative is equal to ^. It is easily 
shown that the condition that the r.h.s. of ( p4[ ) be equal 
to ^ has either two solutions (on the same side of a*j) or 
none. Obtaining these solutions explicitly requires solv- 



ing an algebraic equation of degree six in h - 



Never- 



theless, the conditions for existence of such solutions can 
be written explicitly, as before. In the early shock phase 
there may now be either two, or one or zero fe^'s outside 
of the shock interval. The length of the realizability in- 
terval in the variable u'Q{a^j) is now 0(|^|~^/^). As to 
the coefficients in front of the distributions 5{a — b^), it 
is now found to be 0(|^|~^/^). It follows that 



P 



^'Htt)^\^\- 



^-*±oo. 



(35) 



The time-dependent constant in front the -2 power law 
can again be expressed in terms of the conditional joint 
probability ps^ already introduced, but this is very cum- 
bersome since it involves the solution of the aforemen- 
tioned equation of degree six. 

The theory can be extended to higher-order space 
derivatives but becomes even more cumbersome. Some- 
what superficial inspection (mostly by dimensional anal- 
ysis) indicates that the pdf's have then power-law tails 
with exponent — (3n + 4)/(3n — 1) CU. For even n the 
tail is present for both large negative and positive values. 
For odd n > 2 it is certainly present for large negative 



values and may also be present for positive ones (e.g. for 
n = 3). 

Finally, we turn to the Eulerian time derivative. We 
define 



Pdtuiv^t) = {S (i] - dtu{x,t))) 
From (g) we have 

dtu{x,t) = u'^{a)dta, 



(36) 



(37) 



where dto, is calculated for a given Eulerian position x. 
Time differentiation of a: = a + tu^ (a) gives 



dta - 



uo(a) 



Hence, 



dtu{x,t) = - 



l + tu'^ia)' 
Uo{a)u'Q{a) 



l + tu'aia) 



(38) 



(39) 



Note that the r.h.s. is just —udxU, as could have been de- 
duced from the Eulerian inviscid equation. Substituting 



( |39| ) in ( pq ) and proceeding almost exactly as in Sec. Ill 
we obtain 



pa,u(r?,i)-2\/3t-4£;±(i)|7?r^/^ 



V 



±oo, 



where 



X P4,0 {u, - 



-^,0,u"'\lDa(t) = l) 



(40) 



(41) 



which involves the joint pdf p4.o('«, u' , u" , u"'\1d^u\ = 1) 
of the initial velocity and its first three derivatives at 
an arbitrary location, knowing that this location is not 
within a mature shock interval at time t. Note that the 
pdf of the time derivative (^0|) is just the pdf of the space 
derivative (|3C|) with the change of variable 77 —> — uoC ^^^d 
an extra averaging over uq. This is the result we expect if, 
in (|39|), we neglect the variation of uo{a) near a preshock. 
It is easily shown, when doing the complete asymptotic 
expansion along the same lines as in Sec. 



Ill 



that this 

is indeed the case for the leading-order behavior. This 
theory can again be extended to pdf's of higher-order 
time derivatives which follow the same power laws as for 
space derivatives. 

Obtaining for the pdf of the Eulerian time derivative 
the same law as for the space derivative is not very sur- 
prising. In high-Reynolds number turbulent flows it is 
well known that, when there is a large mean flow, the 
Eulerian temporal structure is, to leading order, deter- 
mined by the spatial structure in the reference frame of 
the mean flow (this is often referred to as the "Taylor 
hypothesis" , but is of course a simple asymptotic re- 
sult). Furthermore, when there is no mean flow, it is 
generally believed that the small-scale temporal struc- 
ture is still determined by the spatial structure, since 



most of the time dependence comes from the sweeping 
of smah-scale eddies by larger energy-containing eddies 
which have much larger but random velocities. For the 
case of burgulence the identical functional forms of ( |30| ) 
and (M) may be seen as a proof of this "random Taylor 
hypothesis" . Note that it is the sweeping by the random 
velocities of the shocks [uo{a) at those locations where 
Mo(a) < 0, UQ(a) = and Mo"(a) > 0] which determines 
the interplay of temporal and spatial structures. Since we 
assumed that the velocity has zero mean value, the ran- 
dom velocities at the shocks can have both signs, so that 
the -7/2 power law tail appears both at large positive 
and large negative values of the Eulerian time deriva- 
tives. Alternatively, one may calculate the pdf of dtu in 
the frame moving with the shock (assuming there is a sin- 
gle shock). In this case one obtains a much steeper law 
C!C \ri\~^. Note that this is not the pdf of the Lagrangian 
time derivative. For unforced burgulence in the inviscid 
limit, this derivative is exactly zero. 



V. VELOCITY INCREMENTS 

We define the Eulerian velocity increment over a sep- 
aration Ax as 

Aue{Ax; X, t) = u{x + Ax, t) - m(x, t). (42) 

Our goal is to find the pdf 

PA«(e, Aa;, t) = {5{e.~ Aue{Ax- x, t))) , (43) 

for values ^ of sign opposite to that of the separation. Let 
us introduce the Lagrangian velocity increment, defined 
outside the Lagrangian shock intervals, as 



Aul{Ax; a, t) — Aue{Ax; a + tuo{a), t). 



(44) 



Proceeding as at the beginning of Sec. Ill, we obtain 

V\u{^,Ax,t) = 



'DL(t) 



S{£_~ AuL{Ax;a,t)) |1 + iiio(a)|da 



(45) 



For given t, Ax and ^, we must now find those Lagrangian 
locations, denoted bk, where the argument of the delta 
function in ( ffSl) vanishes. For this, it is convenient to 
associate to each bk, the point &'j, such that their images 
by the Lagrangian map Ct are separated by a distance 
Aa;, while the velocities differ by ^. We thus have 

uoib'k) ^ uoibk) = ^ (46) 

b'k + tuoib'k) = bk + tuaibk) + Ax. (47) 

The equivalent of (P^ ) is now 



PAu{LAx,t) = 

' {l+tu',{b',)){l+tu',{bk)) 



E 



u'oib',)-u',ibk) 



DlW 



Sia~bk)da). (48) 



We shall here be interested exclusively in situations 
where 



\Ax\ « \t^\ « 1, 



(49) 



which originate from the neighborhood of prcshocks a*j 
where the Taylor expansion (n% may be used |ll| . From 
( p6| ) and (47), we then obtain that bk is a root of the 
following quadratic equation in b : 



(b 



2 


Ax' 


K'(a*j) 



= 0. 



(50) 



We shall see that (pd), together with the realizability 
condition of having real roots not belonging to a shock 
interval, can have either zero, one or two solutions, which 
we shall denote by bjm- We then approximate the pdf 
( [48| ) using the Taylor expansion ( [l2| ) near preshocks, to 
obtain 



P\u{^,Ax,t) 




B, 



5{a — bjm)da 



(51) 



DlW 



where 



[1 + tu'^ja^j) + tu'^'{a^j){bj,n - Q*j)V2]^ 
Bj,ri =t[l + tu'^{a^j) + tuQ {a^j)[bj,ri - a^jf /l] 



A ■ — 



(52) 
(53) 



The realizability conditions associated to (50) lead to 
rather involved expressions for the pdf. Simple scaling 
behavior emerges in two limiting cases. To express the 
corresponding conditions in reasonably compact fashion, 
we shall assume that the third derivative Mo"(a*j) is of 
order unity. 

For \t^\ < |Ax|^/^ the Ajm term in (|l|) dominates. 
For |t^| ^ |Ax|^/'^, the Bjm term dominates. In the 
former case, Ajm can be further approximated 



\A, 



|teP|<'(a.,)(6,™-«*,)l' 



(54) 



(Realizability imposes that ^ and Ax be of opposite sign.) 
Substitution of (|5l) into (|5^) leads essentially to the same 
expression (p^) as for the pdf of velocity gradients, pro- 
vided (i) we replace the Eulerian velocity gradient by 
^/Ax and (ii) multiply the pdf ( PTj) by l/|Aa;|. Hence, 



PAu{^,^X,t) 



^ 



5/2|t|-7/2 



|Ax|^/^|e| 



for |Aa;| < It^l < |Aa;|i/3^ -^ < 0. 



(55) 
(56) 



For 1 > \t^\ > lAxI^/^ the Sj,„ term in (|5^) domi- 
nates. Contrary to the former case, the situation is quite 
different from what has been studied in Sec. |l|. We 
shaU thus give more detaiL With the assumptions made, 
the condition that the quadratic equation (pG) have real 
roots reads 



A{a.j) < -i - <^tr, 



(57) 



which, given the positivity of Ug"(a*j), imphes that t > 
t^j = — l/un(a,f,-). It is then easily checked that one of the 
two roots of (^ is not acceptable because it is within 
the Lagrangian shock interval. The condition that the 
other one be outside this interval reads 



Uq yo-*]) ,2f2 



6 



rr < Wo(a*j)- 



(58) 



Eqs. ( |57| ) and (|5^) play now the role of the realizability 
conditions ( pq ) in Sec. Ill, Proceeding then along the 
same lines as in Sec. Ill, we obtain 



PA«(^,Ax,i)~-F(i)|Aa;||f|, 

for |A.t|1/3 < |<^| < 1, -^ < 0, 
Aa; 



(59) 
(60) 



where 



F{t) ^ / 
Jo 



du"'{u 



lll\2 



P3,0 



--An" 



«Cg(«) 



= 1 



(61) 



and p3,o is defined as in Sec. III. 

We observe that the ^-dependence of our results (|5£ 
and (pfl) for decaying burgulence is essentially the same 
as what was proposed in Ref. ||l| for the forced case. 



VI. CONCLUDING REMARKS 

We have shown here that several results proposed in 
Ref. y for forced burgulence are also valid for decaying 
burgulence and can actually be derived by systematic 
asymptotic expansions, using a Lagrangian approach. 
The results which carry over from the forced to the un- 
forced case are those involving preshocks : -7/2 power law 
for the pdf of velocity gradients and increments and -1-1 
power law for the pdf of increments over suitable ranges. 
For the tail of the pdf of gradients at large positive ^ a 
decaying exponential law of the argument ^^ is generally 
proposed |ll|l^-[l5|]. This result, unrelated to preshocks, 



has no counterpart in the decaying case. Indeed, it fol- 
lows then from D£^/Dt = — C^, where D/Dt denotes the 
Lagrangian derivative, that the pdf of ^ is exactly zero 
for e > 1/t. 

Our results about the -7/2 power law are quite ex- 
plicit : for example we obtain the (time-dependent) con- 
stant D{t), given by (p^), in front of the power law ([3^). 
The expression of D(t) for short times, when mature 
shocks have negligible probability, can be written explic- 
itly in terms of the joint pdf of the first three derivatives 
of the initial velocity at an arbitrary location. For later 
times we need the conditional pdf knowing that no ma- 
ture shock is present at that location. Obtaining this 
exactly, say for Gaussian initial conditions, may be very 
hard. But it is possible to construct lower bounds. For 
example, large deviations theory may be used to show 
that - In D{t) = 0(ln^ t) for i ^ oo |l6|. 

We observe that, formally, our results can easily be 
extended from the case of a space-periodic homogeneous 
initial condition uo{a) (as assumed in Ref. ||l[) to that of 
a random homogeneous mixing initial condition defined 
on the whole real line. For this it suffices to use ergod- 
icity and to replace, in (p^), the mean of the integral 
over the period by limi_^oo(l/(2i)) /_^ . After this, the 
calculation is essentially unchanged. 

It is of interest to point out that without careful han- 
dling of the shock conditions, an incorrect \£,\~^ power 
law is obtained for the left tail of the pdf of gradients. 
Indeed, the pdf (p2) has an overall \^\~^ factor in front 
of the r.h.s. (for large |^|). If in the remaining factor 
we perform the integration over the whole Lagrangian 
interval [0, 1[ without excluding the Lagrangian shock 
interval, we obtain an order unity contribution. It has 
already been pointed out in Ref. Jl^l that a |^|~'^ law 
is obtained from multi-valued solutions of the Riemann 
equation. This is indeed equivalent to using the naive 
Lagrangian map; note that it gives the correct answer 
when using the Zeldovich approximation in cosmology, 
which allows "multi-stream" solutions [[l8|jl9| ]. 

It is shown in Ref. [gO[ that there is a simple relation 
between the velocity gradient and the density of an ad- 
vected passive scalar. When this density is initially uni- 
form, this implies that the power law with exponent —7/2 
also applies to the tail of the density pdf. Analogous re- 
sults can be obtained in several dimensions, where they 
have cosmological implications; this requires the study of 
singularities of multi-dimensional convex hulls p9| . 

Finally, there is a problem which is in a way interme- 
diate between the decaying and the forced case, namely 
"kicked" burgulence. The space-periodic force / appear- 
ing in Im) is then of the form 



/(x,t) = ^/,(x)^(i-t, 



(62) 



where the fj{x) are deterministic or random prescribed 



functions. Between the "kicking times" tj we have decay- 
ing dynamics. At time tj the velocity undergoes a tem- 
poraUy (but not spatially) discontinuous change fj{x). 
When the tj's are equally spaced and all the /j(a;)'s are 
equal, the solution of the Burgers equation converges to 
a space-time-periodic function. Pdf 's obtained by space 
and time averaging have exactly the same scaling prop- 
erties as obtained here for random decaying burgulence. 
Very clean scaling can be obtained by numerically simu- 
lating this problem using a modification of the fast Leg- 
endre transform method of Refs. |p|,pT|. Such questions 
are discussed in Ref. ||^. 
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